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Abstract. We use suspended graphene electromechanical resonators to study the 
variation of resonant frequency as a function of temperature. Measuring the change 
in frequency resulting from a change in tension, from 300 K to 30 K, allows us 
to extract information about the thermal expansion of monolayer graphene as a 
function of temperature, which is critical for strain engineering applications. We find 
that thermal expansion of graphene is negative for all temperatures between 300K 
and 30K. We also study the dispersion, the variation of resonant frequency with 
DC gate voltage, of the electromechanical modes and find considerable tunability 
of resonant frequency, desirable for applications like mass sensing and RF signal 
processing at room temperature. With lowering of temperature, we find that the 
positively dispersing electromechanical modes evolve to negatively dispersing ones. We 
quantitatively explain this crossover and discuss optimal electromechanical properties 
that are desirable for temperature compensated sensors. 
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1. Introduction 

Electronic properties of graphene have been studied extensively [H [2] since the 
first experiments probing quantum Hall effect OH]. In addition to the electronic 
properties, the remarkable mechanical properties of graphene include a high in-plane 
Young's modulus of ~1 TPa probed using nanoindentation of suspended graphene[5], 
force extension measurements [6], and electromechanical resonators [71 El [9]. NEMS 
(nano electromechanical systems) devices using nanostructures like carbon nanotubes 
[TOl [m |12l [131 im [HI US], nanowires [T7] [18] and bulk micromachined structures 
[T9l l20l |2T] offer promise of new applications and allow us to probe fundamental 
properties at the nanoscale. NEMS [22] based devices are ideal platforms to harness 
the unique mechanical properties of graphene. Electromechanical measurements with 
graphene resonators [TUB] suggest that with improvement of quality factor (Q), graphene 
based NEMS devices have the potential to be very sensitive detectors of mass and charge. 
Additionally, the sensitivity of graphene to chemical specific processes [23] offers the 
possibility of integrated mass and chemical detection. The large surface-to-mass ratio 
of graphene offers a distinct advantage over other nanostructures for such applications. 
In order to better understand the potential of graphene based electrically actuated and 
detected resonators [8] and the challenges in realizing strain-engineered graphene devices 
[211 [25] . we experimentally measure the coefficient of thermal expansion of graphene 
{oigrapheneiT)) as a fuuctiou of temperature. Our measurements indicate agraphene(T) < 
for 30 K < T < 300K and larger in magnitude than theoretical prediction [26] . We also 
probe the dispersion, or the tunability, of mechanical modes using the DC gate voltage at 
low temperatures and find that the thermal expansion of graphene, built-in tension and 
added mass play an important role in changing the extent of tunability of the resonators 
[27] and the resonant frequency - parameters that are critical for various applications. 
Additionally, measuring temperature dependent mechanical properties [28] of suspended 
structures down to low temperatures will give insight into strain engineering of graphene 
based devices [241 EH] and also help in understanding the role of rippling in degrading 
carrier mobility at low temperatures. Our experiments probe these issues in detail 
and suggest that monitoring the resonant frequency as a function of temperature could 
provide important information about nanoscale stress that is useful for probing phase 
transition. 

2. Experimental details 
2.1. Device fabrication 

To fabricate monolayer graphene electromechanical resonators, we suspended graphene 
devices using the previously reported [291 (SOl [SH [8] process which starts with 
micromechanical exfoliation of graphene from graphite crystal using an adhesive tape 
on a 300 /im thick degenerately doped silicon wafer coated with 300 nm thick 
thermally grown Si02. Following the location of monolayer graphene flakes using optical 
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microscopy, electron beam lithography is used to pattern resist for fabricating electrodes 
for electrical contact. The electrodes are fabricated by evaporation of 10 nm of chromium 
and 60 nm of gold following patterning of resist. To release the graphene from the Si02 
substrate, a dilute buffered-HF solution is used to selectively etch an area around the 
graphene device by masking the rest of the substrate using polymer resist. Following an 
etch for 5 min 30 sec that results in a 170 nm deep trench in Si02, the device is rinsed 
in Dl water and isopropanol. Critical point drying, to prevent collapse of the device 
due to surface tension, is the final step in fabrication of suspended graphene devices. A 
scanning electron microscope (SEM) image of a suspended graphene device is shown in 
Figure la. 

2.2. Measurement scheme 

The electrical actuation and detection is done by using the suspended graphene device 
as heterodyne- mixer [8l [101132]. The scheme for electrical actuation and detection is also 
shown in Figure la superimposed on the SEM image of the device. We use electrostatic 
interaction between the graphene membrane and the back-gate electrode to actuate the 
motion in a plane perpendicular to the substrate. A radio frequency (RF) signal Vgioj) 
and a DC voltage V^*" are applied at the gate terminal using a bias-tee. Another RF 
signal Vsd{.<^ + is applied to the source electrode Figured^. RF signal applied at 
the gate Vgiu) modulates the gap between graphene and substrate at frequency and 
V^^ alters the overall tension and carrier density in the membrane. The amplitude 
of the current through the graphene membrane at the difference frequency (Aw), also 
called the mixing current Imixi,^^), can be written as [HI Uni [TEl ESI |3l] 

1 

I^U^u) = -^(-^ziuWf'^' + CMuWsnico + Au), (1) 

where G is the conductance of the graphene device, q is the charge induced by the 
gate voltage, Gg is the capacitance between the gate electrode and graphene, z{u) is 
the amplitude of oscillation at the driving frequency u along the ^-axis (perpendicular 
to the substrate). The difference frequency signal (at Au) arises from the product of 
the modulation signals of Vsd{^ + Aw) and G{u). At the mechanical resonance of the 
membrane, the first term in Equation [1] contributes significantly and the second term 
which does not depend on the mechanical motion of the graphene membrane provides 
a smooth background. 

2.3. Results and discussions 

Figure Hb shows the result of such a measurement at 7K for a suspended graphene 
device (Dl) while V^^*" is set at -5V. Using a Lorentzian lineshape for the resonance 
curve, [TOl |35] we can extract the Q ~ 1500 of the resonator (variation of Q with 
temperature is shown in supplementary information). Figure [T}:; shows the colorscale 
plot of /mix (Aw) as a function of V^^*" and / = a;/27r at 7 K. The resonant frequency is 
lowered as the magnitude of is increased - mechanical mode disperses negatively 
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Figure 1. a) Tilted angle scanning electron microscope (SEM) image of a suspended 
monolayer graphene device and the electrical circuit for actuation and detection of the 
mechanical motion of the graphene membrane. The scale bar indicates a length of 
2/^m. b) A plot of the mixing current lmix{^'^) as a function of frequency /(= ^) 
at 7 K with the DC gate voltage V^'~^ = —51/ for device (Dl). The sharp feature in 
the mixing current corresponds to the mechanical resonance, c) Colorscale plot of the 
mixing current as function of frequency / and DC gate voltage V^'~^ at 7 K for device 
(Dl). The dashed line shows the position of the line scan shown in Figure [Tja. The 
dot-dash line indicates the position of the Dirac peak for the graphene device. The 
measured dispersion of this mode of the device is negative. The maxima (blue) and 
minima (red) of the colorbar correspond to 2.6nA and OnA respectively, d) Plot of 
the measured Q for the data shown in Figure [TJ:. The Q drops as IV^^*-^! increases and 
around the Dirac point the Q shows a dip. e) Colorscale plot of Imia;(Aw) as a function 
of frequency / = ^ and the DC gate voltage. V^^ = 13V for device (D2) at 300K. 
Two positively dispersing modes of the device are seen. 
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with IV^^*^!- This is well understood in terms of mode-softening due to the capacitive 
contribution to the energy of resonator [361 118]. The Dirac peak for our device is shifted 
(Vd = 1?>V) from the expected position (0 V) due to unintentional doping during the 
fabrication process 0. Also seen in Figure ^ is the information regarding the amplitude 
of the mixing current. The overall amplitude of the mixing current at resonance scales 
with V^'~' (first term in Equation [1]) and as a result the amplitude of the mixing current 
is very small near V^^ = OV. However, at Vg^^ = 13V, the amplitude of the mixing 
changes abruptly due to the Dirac peak (Vd) pill]. At this point, the term ^ appearing 
in Equation [1] becomes zero. Figure [T]i shows the quality factor (Q) dependance on V^^"" 
calculated from the data shown in Figure [Tfc. The quality factor of the device is largest 
around V^'" = OV and decreases as the gate voltage is swept away from zero. Also, 
at low temperatures, /o decreases with the increase in IV^^*"! and Q follows a similar 
trend. As negative dispersion - caused by softening of spring constant - is accompanied 
by the movement of graphene towards the gate electrode, it leads to an increase in the 
modulated capacitance and results in larger dissipative current [36] . In the neighborhood 
of Vg^'" = Vd, Q is smaller than the value expected from the trend from higher V^'" side. 
One possible reason for this behavior might be that in the vicinity of the Dirac point, 
large device resistance leads to increased ohmic losses [35]. An additional mechanism 
that can result in dissipation is the mechanical deformation modifying the distribution 
of charge puddles of electrons and holes in the graphene sheet. These phenomena 
warrant further detailed investigation [371 [38]. Such a measurement could be used to 
study charge inhomogeneity. Figure [Tfe shows similar measurement for a device (D2) 
at 300K. Here, we observe two mechanical modes which disperse positively because of 
the increase in tension in the graphene membrane due to the electrostatic attraction 
from the back-gate. Most of our graphene devices exhibit an in-built tensile stress at 
300K due to the fabrication process. Multiple ripples on the membrane are also seen in 
Figure la [28]. To better understand the modal dispersion quantitatively, we model the 
graphene membranes within the elastic continuum regime and in a limit where tension 
in the membrane dominates over the flexural rigidity [HI [3, [39]. As a result, the resonant 
frequency /o of the graphene membrane can be written as 



where L is the length of the membrane, w is the width, t is the thickness, p is the 
mass density, To{T) is the in-built tension and F is the tension at a given temperature 
T and V^^"". The functional form of r(ro(T), V^^*") is dependent on the details of the 
model used to take into account the electrostatic and elastic energies. (Our model is 
described in detail in the supplementary information.) Resonant frequency at zero gate 

voltage is given by /o(0) = '■, therefore, from /o(0) an independent estimation of 



2L y ptw 

f A shift in the Dirac peak away from Vj^'~^ ~ OV is desirable in our experiments because the actuation 
efficiency is feeble at Vj-"^ = OV and that would make the observation of physics near the Dirac point 
inaccessible for electromechanical measurements. 
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ro(T) and p is not possible. However, resonant frequency dispersion with gate voltage 
foiy^^) allows us to estimate p and Tq{T) simultaneously. Using fits based on the 
continuum model (calculations described in supplementary information) we estimated 
the mass density p = T Apgraphene and the in-built tension ro(300i^) to be 68.6 nN for 
the fundamental mode, where Pgraphene is the mass density of pristine graphene for the 
device data shown in Figure [1^. We attribute this extra mass and in-built tension to the 
resist residue that can get deposited on the graphene membrane during the fabrication 
process and these values of mass density are similar to ones reported by Chen et al. [8]. 
We note that the dispersion of the two modes shown in the Figure [1^ is different and 
our modal calculations give p = 7.Spgraphene and ro(300-ft') = 107.6 nN for the upper 
mode. However, the simple model of a rectangular membrane predicts higher order 
modes at much higher frequencies - this is a limitation of the simple assumptions we 
have made. The presence of other resonant modes at small intervals of frequency / seen 
in Figure [1^ (also seen in Figure H^) indicates a deviation from this simplified picture 
of membrane under tension due to either edge-modes [9], or due to the rippling/curling 
of the membrane [28l HQ]. In the later part of this report, we discuss in detail the 
reasons for positive or negatively dispersing modes and temperature induced change in 
dispersion. 

3. Probing thermal expansion of graphene 

We now consider how the resonant frequency (/o) evolves as a function of the 
temperature. Figure shows the result of an evolution of a mode as a function 
of temperature at V^^*" = 15 V The resonant frequency increases as the device is 
cooled from room temperature. This increase has been seen in all the devices we have 
studied. The degree of frequency shift varies from one device to another depending on 
the device geometry. The origin of this frequency shift with temperature is the increase 
in tension in graphene due to the expansion/contraction of substrate, gold electrodes and 
graphene. The frequency shift can be understood by taking into account the contribution 
of various strains as the device is cooled below 300 K. The three main contributions are 
- firstly, the thermal strain in unconstrained graphene egraphene(T) = Jj, agrapheneit)dt 
[§| due to the coefficient of thermal expansion of graphene agraphene{T), second, the 
thermal strain due to the gold electrodes egoidiT) = fj, agoid{t)dt and lastly the 
contribution of the strain induced by the substrate egubstrateiT)] here agoidiT) is the 
coefficient of thermal expansion for gold [41J. The strain in gold electrodes plays 
an important role due to the geometry of the device. The under-etch that releases 
the graphene membrane also etches under the graphene covered by the electrodes - 
resulting in the graphene membrane being suspended off the gold electrodes [8] (see 
supplementary information for SEM image). The geometry of the resulting device 

X Data acquisition was done during a single sweep over twelve hours to allow the resonator to equilibrate 

and the window of acquisition window automatically adjusted to follow the resonance. 

§ All through this report we measure the strain relative to the strain at room temperature (300 K). 
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Figure 2. a) Plot showing tlie evolution of the resonant frequency of a mode for device 
(D2) as a function of temperature for V^'-^ = 15 V. Inset shows the schematic of all 
the strains external to the suspended graphene membrane as the device is cooled below 
300 K. b) The plot of expansion coefhcient of graphene as a function of temperature. 
Data from two different devices together with theoretical prediction of N. Mounet et 
al. [2^. The shaded area represents the errors estimated from the uncertainty of the 
length of the flake, width of the electrode and Young's modulus of graphene. 

is shown in the inset to Figure [2^. The elastic strain in the curved substrate can 
be calculated using Stoney's equation |42]. Its contribution to the net strain in the 
graphene membrane is very small (see supplementary information) and therefore upon 
cooling, change of tension in the graphene is due to contraction of gold electrode and 
expansion/contraction of graphene. We assume that Young's modulus of graphene does 
not vary significantly over the temperature range |13] . At the interface of gold electrodes 
supporting the graphene membrane, the net force must balance to zero; however, the 
stresses are different considering the crossectional area of gold electrodes (~500nm x 
60nm) and graphene (~500nm x 0.3nm). The large difference in the cross-sectional area 
implies that the effective stiffness of gold electrodes is large compared to the stiffness of 
graphene. As a result, to a very good approximation, the total elastic strain at a given 
temperature in graphene that is confined by "rigid" gold electrodes is egrapheneciamped = 
(^graphene{T) + esubstrate(T) - egoid{T)^^^^i^f^ , where Weiectrode IS the average of the width 
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of gold electrodes holding the suspended flake. The change in tension in the membrane 
can be written as a function of temperature as Aro(T) = wtegraphenedamped{T)E graphene, 
where Egraphene is Young's modulus of graphene. Measuring the tension ro(T) as a 
function of temperature offers a way to track the thermal strain in graphene membrane. 
Figure [2^ shows the evolution of resonant frequency as a function of temperature from 
a device (D2) at V^'" = 15 V where the increase in frequency is largely due to the 
contraction of the gold electrodes. However, this rate of increase of resonant frequency 
is significantly reduced due to the negative agraphene for all T< 300 K from the case 
of frequency change including only gold's contraction. Using such a measurement 
of frequency shift while assuming uniform expansion of all the materials and using 
Equation [21 it is possible to calculate the expansion coefficient from the frequency at 
Vg^^ = V as 

P d , /rri\ /rri\^ electrode \ /Q\ 

CXgrapheneyJ- ) = — ^JolUJ X — h — (^Csufesfrate U ) — ^gold\-^ ) 7 ) W) 

W-t -t^ graphene 1j 

Figure [2]d shows the result of calculating agraphene for two devices using this analysis 
and comparison with the theoretical calculation for agraphene by N. Mounet et al. 
. We find that agraphene is negative and its magnitude decreases with temperature 
for T< 300 K. At room temperature, agraphene ~ —7 x 10~^K~^ , which is similar 
to the previously reported values measured by others [H [28]. At 30 K, agraphene ~ 
— 1 X lO^^i^^^. The deviation of agraphene from the theoretically predicted values 
can possibly be due to the presence of the impurities on graphene membrane. The 
knowledge of agraphene is essential for the fabrication of the devices intended for 
strain engineering applications [211 EH]. Strain engineering of graphene devices at 
low-temperatures using this picture can improve device performance, for example by 
enabling temperature compensation [44]. A simple design rule for width of electrodes to 
achieve temperature compensation in the vicinity of a temperature Tq is to ensure that 
the w electrode = L X I o'araphenciTo) i ^^gjj^g equation [3]) . Additionally, these measurements 
indicate that using NEMS resonators for measuring internal stress of nanowires of metal 
and phase change materials as a function of temperature, using phase locked loop (PLL) 
technique [21], can provide useful information about stresses and strains in individual 
nanostructures. 

Here, it is important to point out that, residue on the graphene sheet can also affect 
the calculations of the coefficient of thermal expansion. This can possibly be one reason 
for its deviation from the theoretically predicted values. In our calculation, we have 
assumed that presence of residue does not affect the expansion of graphene, which may 
not be the complete description of the experimental system. However, these values of 
agraphene stiU remain valuable for device engineering as the complete elimination of the 
resist residue from the device is difficult. 
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Figure 3. a) Measured dispersion of an electromechanical mode as function of 
temperature for device (D3). b) Closeup of the experimentally measured low- 
temperature data showing the non-monotonic dispersion. c) Modeling of the 
dispersion, incorporating thermal expansion of graphene, at various temperatures 
by varying temperature and the parameter A. d) Closeup of the calculated low- 
temperature data showing the non-monotonic dispersion similar to the experimentally 
measured dispersion. 



4. Modal dispersion of graphene NEM resonator 

We next consider the dispersion of the resonant modes with temperature. Figure [3^ 
shows the dispersion of resonant mode for a device(D3) at various temperatures down to 
6K. Upon cooling, the resonant frequency increases and its tunability using gate voltage 
is reduced, an undesirable feature. Further, at low temperatures, resonance frequency 
first decreases with IV^^*"! and then it increases (Similar evolution of modal dispersion 
with V^^ has been observed by Chen et al. [8].). We now try to understand the 
evolution of the modal dispersion as a function of temperature. The two limiting cases 
of modal dispersion - pure positively and negatively dispersing modes can be respectively 
understood by considering the limits where the tension (or in the case of flexural modes 
mechanical rigidity) dominates with increasing V^^ and the case when capacitive pulling 
leading to the softening of the spring constant dominates [361 [HI [35] . A simple model [18] 
to understand the intermediate regime where these two interactions compete is to model 
the resonator with an intrinsic spring constant Ki = k+a{V^'")'^+P{Vj^'")'^+H.O.{Vj^'") 
with a and /3 as constants of the system (when /3 > 0, this is consistent with the 
positive dispersion of modes with an increased V^'-^). A second contribution due to the 



Graphene NEMS resonator 



10 



electrostatic interaction (capacitive coupling), softens the intrinsic spring constant Ki 
so that the effective spring constant 

K.u - K. - l(Vrf^. (4) 

with Cg being the capacitive coupling of the membrane with gate electrode. As V^'" is 
varied the K^ff varies and the modal dispersion changes from a negatively dispersing 
mode to a positively dispersing mode at large Vg^'-^ - the value of V^*-^ at which the 
crossover happens is a measure of the relative contribution of capacitive and elastic 
energies. This phenomenological model describes the generalized modal dispersion 
and can be connected to the properties of the resonator (described in supplementary 
information), a and /3 depend on temperature through the coefficient of thermal 
expansion, and also on the elastic constant of graphene; as a result they change with 
temperature, leading to the evolution of modal dispersion as a function of temperature. 
Additionally the capacitive contribution to energy is expected to change as a function 
of temperature. The observed modal dispersion behavior at all temperature can be 
quahtatively understood from the contraction of the suspended gold electrodes with 
temperature, which leads to increased tension in the membrane. This increase in in- 
built tension is accompanied by the reduction in tunability with V^'-' (model described 
in the the supplementary information). Results of our calculation (Figure [St) clearly 
show that we can successfully model the temperature evolution of dispersion, using 
only a single fit parameter in our model. It is critical to understand the origin of 
observed dispersion of resonant frequency as a function of temperature due to the 
desirable property of resonator - its tunability. Additionally, if the loss mechanisms 
are frequency independent such tunability can increase the Q of the system |35j . 

In order to use graphene electromechanical resonators for applications like mass 
and charge sensing it is important to understand the microscopic origin of the modes 
and to test the suitability of continuum models [3, [8], [39] beyond basic properties such as 
dispersion. Scanned probe measurements to image modes [9] indicate that the origin of 
these modes is complex - for instance some of the modes are associated with the edges 
of the graphene membrane. The reason this is critical is that the notion of effective mass 
of mechanical modes [221 US] is associated with the spatial distribution of the mode [jj]. 
We try to understand the nature of modes in graphene resonators by studying devices 
with multiple resonances in the measurement frequency window. Figure [1^ shows two 
higher order modes along with fundamental mode for a graphene resonator at 7 K. 
The dispersions of these three modes are different. Though it is difficult to determine 
the exact nature of these modes using our technique, a semi-quantitative explanation is 
provided here. At lower temperatures, negative dispersion can be understood in terms of 
spring constant softening as discussed earlier. Following Equation [U the effective spring 

II The effective mass meff of a mode is dependent on the mode in flexural oscillations. For the case of 
string under tension the me/ / is independent of the mode [3S] . 
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constant -ft"e// can be used to calculate the modal dispersion with gate voltage (V^^*"), 



where /o = ^J:;^^ and me// is the effective mass for a given mode 



Using Equation |5l we can describe the higher modes also by using /o and modal 
mass {nieff) as the fitting parameters. For a sheet under tension with uniform loading, 
the effective modal mass (^e//) for different modes does not change, unlike the case 
of flexural modes, and is equal to riieff = 0.785pLwt |15]. However, our fitting gives 
different effective masses for different modes, which is expected as these are not the 
higher order harmonics of the fundamental mode (in the model described above, the 
higher order modes will be integral multiples of the fundamental mode). This suggests 
that the simple picture of rectangular membrane under tension within the continuum 
description does not work well to describe the system. This could be due to four main 
reasons: a) non-uniform loading of the membrane (due to resist residue) [8] can modify 
the meff for different modes, b) due to the presence of edge dependent modes |9|, c) 
due to the curvature of membrane or rippling of graphene (as seen in Figure 
[IK) the effective stiffness of the modes could be a very complex quantity with a tensor 
nature and could be significantly different from the ideal value of ~ 1 TPa as observed 
for the case of rippled carbon nanotubes |l6l H?! HH] and d) the capacitance to gate 
(Cg) for the graphene membrane is likely to be mode dependent. Further experiments 
with pristine unrippled graphene resonators are needed to clarify our understanding of 
the microscopic origin of modes. The presence of multiple closely spaced modes can 
be potentially useful for mass spectrometry as the position of the added mass can be 
extracted accurately 



5. Conclusions 



To conclude, our measurements of Q, for graphene resonators, as a function of DC gate 
voltage indicate that larger dissipation in the resonator can occur around the Dirac 
point. Using these NEMS devices measurement of coefficient of thermal expansion of 
graphene at low temperatures is possible. Our measurements indicate that agraphene 
is negative for all temperatures between 300 K and 30 K and larger in magnitude 
than the numbers predicted by theoretical calculations [26]; this could be critical in 
designing strain engineered devices using graphene [211 125]. Additionally, the modal 
dispersion of graphene resonators is affected by the thermal expansion of graphene 
reducing the tunability. The continuum description of the mechanics of rectangular 
graphene membrane is inadequate for explaining the resonances due to the presence of 
non-uniform impurities and rippling of the membrane. 
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Figure 4. a) Colorscale plot oilmix{^^) showing tliree negatively dispersing modes 
with varying slopes, b) Fitted data for varying effective mass for the three modes. The 
inset shows two possible modes for a uniform rectangular membrane under tension. The 
effective mass is expected to be independent of the mode - an observation that is not 
experimentally seen within the measured frequency range. 
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